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Abstract 

In the framework of the coupled channel model the mass shifts of the P-wave excitations 
of Ds and Bg mesons have been calculated. The corresponding coupling to DK and BK 
channels is provided by the effective chiral Lagrangian which is deduced from QCD and does 
not contain fitting parameters. The strong mass shifts down for 0"'' and 1"*" states have been 
obtained, while 1+ and 2+ states remain almost at rest. Two factors are essential for large 
mass shifts: strong coupling of the 0"*" and 1+ states to the S-wave decay channel, containing 
a Nambu-Goldstone meson, and the chiral flip transitions due to the bispinor structure of 
both heavy-light mesons. The masses M{B*{0+)) = 5710(15) MeV and Af(B^(l+')) = 
5730(15) MeV are predicted. Experimental limit on the width r{Dsi(25S6)) < 2.3 MeV 
puts strong restrictions on admittable mixing angle between the 1+ and 1+ states. 

1 Introduction 

The heavy-light (HL) mesons play a special role in hadron spectroscopy. First of all, a HL 
meson is the simplest system, containing one light quark in the field of almost static heavy 
antiquark, and that allows to study quark (meson) chiral properties. The discovery of the 
Ds{2317) and Z?<j(2460) mesons [D 12] with surprisingly small widths and low masses has given 
an important impetus to study chiral dynamics and raised the question why their masses are 
considerably lower than expected values in single channel potential models. The question was 
studied in different approaches: in relativistic quark model calculations [3]-[6], on the lattice 
[7j, in QCD Sum Rules [l [9], in chiral models [IO]-[l2] (for reviews see also [131111]). The 
masses of Ds{0~^) and Ds(l~^ ) in closed-channel approximation typically exceed by ~ 140 and 
90 MeV their experimental numbers. 

Thus main theoretical goal is to understand dynamical mechanism responsible for such large 
mass shifts of the O"*" and l"*" levels (both states have the light quark orbital angular momentum 
I = and j = 1/2) and explain why the position of other two levels (with j = 3/2) remains 
practically unchanged. The importance of second fact has been underlined by S.Godfrey in [5]. 

The mass shifts of the Ds{0'^,l^ ) mesons have already been considered in a number of 
papers with the use of unitarized coupled-channel model [TS], in nonrelativistic Cornell model 
|16j . and in different chiral models [17j-[19j. Here we address again this problem with the aim 
to calculate also the mass shifts of the Ds{l~^ ) and Bs{0~^ , 1"*" ) states and the widths of the 

and 1"*" states, following the approach developed in [18], for which strong coupling to the S- 
wave decay channel, containing a pseudoscalar (P) Nambu-Goldstone (NG) meson, is crucially 
important. Therefore in this approach principal difference exists between vector-vector (1^1^) 
and VP (or PP) channels. This analysis of two-channel system is performed with the use of 
the chiral quark-pion Lagrangian which has been derived directly from the QCD Lagrangian 
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|20j and does not contain fitting parameters, so that the shift of the D*{0'^) state ~ 140 MeV 
is only determined by the conventional decay constant fx- 

Here the term "chiral dynamics" implies the mechanism by which in the transition from one 
HL meson to another the octet of the NG mesons (p is emitted. The corresponding Lagrangian 

^Lfcm = q{(yr) exp(i75(/>//^)g, (1) 

contains the light-quark part, exp(i75(/)//^), where (j) is the SU{2>) octet of NG mesons and the 
important factor 75 is present. In the lowest order in (j) this Lagrangian coincides with well- 
known effective Lagrangian ALgfj suggested in |21j.[22j. where, however, an arbitrary constant 
QA is introduced . At large Nc, as argued in [21], this constant has to be equal unity, qa = ^■ 
In [101 [T71 [22] this effective Lagrangian was applied to describe decays of HL mesons taking 
QA < 0.80. 

More general Lagrangian ALpcM ^ was derived in the framework of the field correlator 
method (FCM) [201 [23], in which the constant = 1 in all cases, and which contains NG 
mesons to all orders, as seen from its explicit expression ([T|). 

In [21] with the use of the Dirac equation it was shown that in the lowest order in (p 
ALpcM = A-^eff) if indeed gA = ^- In our calculations the ALpcM was used to derive the 
nonlinear equation for the energy shift and width, AE = AE — as in |T8]. We do not assume 
any chiral dynamics for the unperturbed levels, which are calculated here with the use of the 
QCD string Hamiltonian [251 [26], because the mass shift AE appears to be weakly dependent 
on the position of unperturbed level. Nevertheless, the uncertainty in the final mass values 
is due to a poor knowledge of the fine structure (FS) interaction in the initial (unperturbed) 
P-wave masses. 

It is essential that resulting shifts of the J^(0^,1^ ) levels are large only for the Ds,Bs 
mesons, which lie close to the DK, D* K, BK, B* K thresholds, but not for the D{1P), B{1P) 
mesons, in this way violating symmetry between them (this symmetry is possible in close- 
channel approximation) . In our calculations shifted masses of the -Ds(0+) and 5^(0+) practically 
coincide with those for the D{0~^) and -6(0"'"), in agreement with the experimental fact that 
Mc^p{D{0+)) = 2350 ± 50 MeV [27] is equal or even larger than Mexp(Ds(0+)) = 2317 MeV. 
The states with j = 3/2 Ds{l~^, 2"*") and -0(1"*", 2"*") have no mass shifts and for them the mass 
difference is ~ 100 MeV, that just corresponds to the mass difference between the s and light 
quark dynamical masses. 

For the £'^(1"*' ) and 5^(1+ ) mesons calculated masses are also close to those of the D and 
B mesons. Therefore for given chiral dynamics the J^(0"^, 1"^ ) states cannot be considered as 
the chiral partners of the ground-state multiplet J^(0~, 1^), as suggested in [llj. 

We also analyse why two other members of the IP multiplet, with = 2+ and 1+, do 
not acquire the mass shifts due to decay channel coupling (DCC) and have small widths. Such 
situation occurs if the states 1"*" and 1"*" appear to be almost pure j = | and j = ^ states. 
Still small mixing angle between them, \(j)\ < 6°, is shown to be compatible with experimental 
restriction on the width of Dsi(2536), admitting possible admixture of other component in the 
wave function (w.f.) < 10%. 

In our analysis the 4-component (Dirac) structure of the light quark w.f. is crucially im- 
portant. Specifically, the emission of a NG meson is accompanied with the 75 factor which 
permutes higher and lower components of the Dirac bispinors. For the j = 1/2, P -wave and 
the j = 1/2, S -wave states it is exactly the case that this "permuted overlap" of the w.f. is 
maximal because the lower component of the first state is similar to the higher component of 
the second state and vice versa. We do not know other examples of such a "fine tuning". On 
the other hand in the first approximation we neglect an interaction between two mesons in the 
continuum, like DK, etc. 

In present paper we concentrate on the P-wave B, Bg mesons and the effects of the channel 
coupling. While the IP levels of the D,Ds mesons are now established with good accuracy 
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[I] ) [2] ) |2Tj , for the B, Bg mesons only relatively narrow 2"'", l"*" states have been recently observed 
[28], [29]. According to these data the splitting between the 2"*" and l"*" levels is small, ~ 20 — 10 
MeV, while the mass difference between Bs{2'^) and B{2'^) states is large ~ 100 MeV, as for 
the Ds{2+) and D{2+) mesons. 

The actual position of the B{lP),Bs{lP) levels is important for several reasons. Firstly, 
since dynamics of [qb) mesons is very similar to that of qc^ the observation of predicted large 
mass shifts of the -65(0"^, 1"^ ) levels would give a strong argument in favour of the decay channel 
mechanism suggested here and in [TS]. It has been shown in [3^ that the mass of Bs{Q'^) can 
change by 150 MeV in different chiral models. Secondly, experimental observation of all P- 
wave states for the B,Bs mesons could clarify many unclear features of spin-orbit and tensor 
interactions in mesons. Understanding of the decay channel coupling (DCC) mass shifts could 
become an important step in constructing chiral theory of strong decays with emission of one 
or several NG particles. 



2 Mixing of the 1+ and 1+ states 

It is well known that in single-channel approximation, due to spin-orbit and tensor interactions 
the P-wave multiplet of a HL meson is split into four levels with = 0+, i;^, 1^, 2+ [ST]. Here 
we use the notation H(L) for the higher (lower) l"*" eigenstate of the mixing matrix because a 
priori one cannot say which of them mostly consists of the light quark j = 1/2 contribution. 
For a HL meson, strongly coupled to a nearby decay channel (DC), some member (s) of the 
P-wave multiplet can be shifted down while another not. Just such situation takes place for 
the -Ds(lP) multiplet. The position of the levels with J = |, which remains unshifted, will be 
important in our analysis. 

The scheme of classification, adapted to a HL meson, in the first approximation treats the 
heavy quark as a static one and therefore the Dirac equation can be used to define the light 
quark levels and wave functions [lOj. Starting with the Dirac's P-wave levels, one has the states 
with j = 1/2 and j = 3/2. Since the light quark momentum j and the quantum number x are 
conserveqj, they run along the following possible values: 
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The HL meson w.f. can be expressed in terms of the light quark w.f. 

.m. 



the Dirac bispinors 



|cT) + C 



1.1 

" 2 ' 2 ' 



where C/iMi j2M2 corresponding Clebsch-Gordan coefficients. 



(3) 
(4) 



Later in the w.f. we neglect possible (very small) mixing between the D{1-^^^), D{1^^^) states 



and also between Ds{2'^^^), Ds{2'^^^) states; however, physical Ds{l^) states can be mixed via 



^we use here the standard notation >f = Tli + |l for ~ |^ 



3 



open channels and tensor interaction. The eigenstates, defining the higher IJj and lower 1^ 
levels, can be parameterized by introducing the mixing angle (p: 

1 3 

\'^h) = cos (t)\j = -)+ sin ct>\j = -), (5) 

and 

1 3 

|1+) = -sin(/)|j = -) + cos(/)|j = -). (6) 

Later we will show that just the level with small |0| ^6°, being almost pure j = | state, has 
no DC (hadronic) mass shift. In opposite case when 1^ is mostly J = f state it is convenient to 
redefine in the equations @, ^ the mixing angle as (/> — > 90° — (p, performing similar analysis. 

In general, the structure of the mixing is important because it defines the order of levels, 
the mass shift down of the l"*" state, as well as the mass shift and the width of another 1+ 
level. One of our goals here is to understand why if the coupling to nearby continuum channel 
is taken into account, the position of the 2+ and 1^ levels does not change (within 1-3 MeV) 
while 0^, 1^ levels acquire large DC shifts. 



3 Chiral Transitions 



To obtain the mass shift due to DCC effect we use here the chiral Lagrangian ([T]), which includes 
both effects of confinement (embodied in the string tension) and Chiral Symmetry Breaking 
(CSB) (in Euclidean notations): 



LFCM = i / d^xip'^id + m + M)^ 



(7) 



with the mass operator M given as a product of the scalar function W{r) and the SU(3) flavor 
octet, 



M = W{r) exp(z7; 



where 



/ 4 



ipaXa = V2 



V2 

K- 



+ 



2ll 



K+ 



(8) 



(9) 



Ve / 



Taking the meson emission to the lowest order, one obtains the quark-pion Lagrangian in the 
form 

ALfcm = - j ^j+(x)(T|x|75^^^-(/'fc(a;)d*^3;. (10) 

Writing the equation as ALpcM = — J Vifdt, one obtains the operator matrix element 
for the transition from the light quark state i (i.e. the initial state i of a HL meson) to the 
continuum state / with the emission of a NG meson {(pa^a)- Thus we are now able to write 
the coupled channel equations, connecting any state of a HL meson to a decay channel which 
contains another HL meson plus a NG meson. 

In the case, when interaction in each channel and also in the transition operator is time- 
independent, one can write following system of equations (see [32] for a review) 



[{Hi-E)5u + Vu]Gif = l. 



(11) 



Such two-channel system of the equations can be reduced to one equation with additional DCC 
potential, or the Feshbach potential [33], 



(Hi - E)Gu - Vi 



1 



12 



Ho — E 



:V2iGu 



1. 



(12) 
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Considering a complete set of the states |/) in the decay channel 2 and the set of unperturbed 
states \i) in channel 1, one arrives at the nonlinear equation for the shifted mass E, 

E = eP - Y,{i\V,2\f)^ (/|^2i|i). (13) 

Here the unperturbed values of E^f^ are assumed to be known beforehand, while the interaction 
Uif is defined in (llOp . A solution of the nonlinear equation (jlSh yields (in general a complex 
number E = E — ^) one or more roots on all Riemann sheets of the complex mass plane. 



4 Calculation of the DCC shifts 

To calculate explicitly the mass shifts, we will use the Eq. (jl3p in the following form: 



= - V '^J'^'^?:/ , (14) 

^ Ef- mh] 



where m('^)[z] is the initial mass, m[i] - is the final one, Ef = cvd + u)k is the energy of the final 
state, and the operator V provides the transitions between the channels (see the comment after 
Eq. (Ho])). 

In our approximation we do not take into account the final state interaction in the DK 
system and neglect the D-meson motion, so the w.f. of the i, /-states are: 

\f >=-^k{v)®-^d{Mj), \i>=-^DAMi), (15) 

where 

gipr 

*'^<'" - TSmS '''' 

is the plane wave describing the X-meson and ^oi^f)-, ^oA^i) ^^^^ the HL meson w.f. at 
rest with the spin projections My, Mj, respectively. 
We introduce the following notations: 



^k{v>) = + W£,(p) = Y'p2 + (17) 

so that in the final state the total energy is Ef = uj£) + ujk, while 

Tf = Ef — rriD — rriK (18) 

is the kinetic energy. Also it is convenient to define other masses with respect to nearby 
threshold: mthr = + mo, 

Eq = m^^^[Ds]—m£)—mK, 6m = m[Ds] — m^'^^[Ds], IS. = Eo+Sm = m\Ds\—mD—mK, (19) 

where A determines the deviation of the Dg meson mass from the threshold. In what follows 
we consider unperturbed masses niQ ( J^) of the {Qq) levels as given (our results do not change 
if we slightly vary their position, in this way the analysis is actually model-independent). 
Using these notations, the Eq. (ll3p can be rewritten as 

- A = .F(A), (20) 

where 

d3p ^\{m, V p.M,]" 



' (2vr)^ ^ Tf{p) - A 
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Mi 



V 



and 

■Mf) = - hl>^oSM,)a\v\^,y-^D{Mj)^===d'v, (22) 

The function ^(A) for negative A diminishes monotonously so there exists a final (critical) 
point, 

= ^(-0). (23) 

Thus, while solving the Eq. fpU]) . one has two possible situations: Eq < Eq^^^ and £"0 > E^"^ 
(see Fig. [IJ. 

In the first case Eq. (j20p has a negative real root A < and the resulting mass of the 
Dg meson appears to be under the threshold. In the second case Eq. ()20p has a complex root 
A = A' + zA" with positive real part A' > and negative imaginary part A" < 0. To find latter 
solutions one should make analytic continuation of the solution(s) from the upper halfplane of A 
under the cut, which starts at the threshold, to the lower halfplane (second sheet). This solution 
can be also obtained by deforming the integration contour in Tf{p). In actual calculations we 
take infinitesimal imaginary part A", proving that A does not change much for finite A" (the 
similar procedure has been used in [H]). Finally, the resulting mass of the Dg meson proves to 
be in the complex plane at the position A' — i|A"|, i.e. the meson has the finite width F = 2A". 

For further calculations we should insert the explicit meson w.f. to the matrix element (|22|) . 
As discussed above, in a HL meson we consider a light quark q moving in the static field of a 
heavy antiquark Q, and therefore its w.f. can be taken as the Dirac bispinor: 



00 



v^,'"=(, ,,£r:';' . /(/=+s=).-=*=i, (24) 

\^(-l) 2 f{r)VLjVM) ^ 
where the functions g{r) and /(r) are the solutions of the Dirac equation: 

g' + -{Eg + mq + U -Vc)f = 0, 

/' + f + {E,-mg-U -Vc)g = ^. 



(25) 



Here the interaction between the quark and the antiquark is described by a sum of linear scalar 
potential and the vector Coulomb potential with = const: 

U = ar, Vc = --, P=^as. (26) 

r 6 

Introducing new dimensionless variables 

x = r^/a, eg = Eg/^/a, ^ig = mq/^/a, (27) 
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Figure 2 



and new dimensionless functions 
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a 



a 



3/4 



F{x) 



oo 

J {F^ + G^) = 1, 



(28) 



we come to the following system of equations: 



G' + -G- 

X 

F'-^F + 

X 
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e„ + + X -\ — 

X 

_ _ P 

En Un X -\- 

X 



F = 0, 
G = 0. 



(29) 



This system has been solved numerically. 
Using the parameters from the papers 

= 0.18 GeV^ as 
uis = 210 MeV, ruq 

we obtain the following Dirac eigenvalues e: 



= 0.39, 
4 MeV, 



(30) 





Qq, fig = 0.01 


Qs, fis = 0.5 


-1 


1.0026 


1.28944 


+1 


1.7829 


2.08607 


-2 


1.7545 


2.08475 



(31) 



and corresponding eigenfunctions G, F are given in Fig. O 

Our choice of a and as is a common one in the frame of the FCM approach, and the 
value of the light quark mass really does not influence here on any physical results because 
of its smallness in comparison with the natural mass scale ^/o^. The strange quark mass is 
taken from [35], where it was found from the ratio of experimentally measured decay constants 
f{Ds)/f{D); the same value can be obtained by a renormalization group evolution starting 
from the conventional value ms{2 GeV) = 90 it 15 GeV. 

Later we use the simplified notations for the quark bispinors: 



Now, using explicit expressions for the spherical spinors, 



def 



(32) 



l+l/2,l,M 



i±M.v 



-2j~ ^/,M+l/2 



l~l/2,l,M 



2j+2 ^l,M-l/2 



j±M+lv 
2j+2 J^/,M+l/2 



(33) 
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Figure 3: $0,2(9) functions 
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and the expansion : 



(34) 



after cumbersome transformations (which are omitted in the text) we obtain the transition 
matrix elements: 



V12 



Mi,M2 



V'UMi) c7|r|75^ V'2(M2) — === 



(35) 



V32 



■ / V';^(M3)(j|r|75^V'2(M2 

J IK 



M3,M2 



VVrKp) 




-1 +2 " 

-\/2 +^/3 

-\/3 +^/2 

-2 +1 



where 



$0(9)= / jo{qx)xdx Gi{x)F2{x) - Fi{x)G2{x) 



00 

!(Q') = y 32{qx)xdx 



G3(x)F2(x)-F3(x)G2(x) 



(36) 



(37) 



Notice that because of different signs of the -Fi(x) and ^2,3(2;) functions (while the (^1^2, 3 func- 
tions are all positive) on almost all real axis, the integral ^2 appears to be strongly suppressed 
in comparison with the integral <I>o- This fact is confirmed by numerical simulations (see Fig. 

ED. 

Finally, introducing universal functions 



-^0,2 (A) 



a f p{Tf)iOD{Tf)dTf y 



2-71^ fl J Tf + niD + mK Tj - A 



(38) 
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Table 1: ^^(O"'') -meson mass shift due to the DK decay channel and i?s(0"*')-meson mass shift 
due to the BK decay channel (all in MeV) 



state 




j^(theor) 


^(oxp) 


6m 




2475 (30) 


2330(20) 


2317 


-145 


i?.(0+) 


5814(15) 


5709 (15) 


not seen 


-105 



we come to the following equations to determine meson masses and widths: 



Ds{lt) 



ii;o[0+]-A = ^o(A), 

E^[l+] - A = cos2 . ^o(A) + sin2 • ^2 (A), 
Eq[1-^] -A' = sin2 (j) . ^o(A') + cos2 • ^2 (A') 
r[l+] = sin2 (j) . fo(A') + cos2 cj) . f2(A'), 
ii;o[2+2]-A' = ^-^2(A'), 

n2y = l-r,{A'). 



(39) 



5 Results and discussion 

In this chapter, using the expressions ([39]) to define the Dg and Bg meson mass shifts, we 
present and discuss our results. We will take into account the following pairs of mesons in 
coupled channels (i refers to first (initial) channel, while / refers to second (decay) one): 



i 




^s(l+) 




f 


D{0-)+K{0') 


D*{1^)+K{0-) 


D*{l-)+K{{)-) 


i 


Bs{0+) 




Bs{2+) 


f 


B{0-) + K{0-) 


B*{l-)+K{0-) 


B*{l-)+K{{)-) 



(40) 



In our calculations we use the following meson masses and thresholds (in MeV): 



mj^+ = 1869, 
mD*+ = 2010, 
m^+ = 5279, 
ms* = 5325, 



m£)+ -|- mj^- = 2363, 
m£).+ -|- m^- = 2504, 
mB+ + mj^- = 5772, 
ms* + mrix- = 5819. 



(41) 



The results of our calculations are presented in Tables dHSl A priori one cannot say whether 
the |j = ^) and \j = |) states are mixed or not. If there is no mixing at all, in this case the 



width r(Z)si(2536)) = 0.3 MeV is obtained in [36], while the experimental limit is T < 2.3 MeV 
|27] and recently in [37] the width T = 1.0 ± 0.17 MeV has been measured. Therefore small 
mixing is not excluded and here we take the mixing angle (j) slightly deviated from = 0° (no 
mixing case). Then we define those angles (j) which are compatible with experimental data for 
the masses and widths of both 1"*" states. 

The small value (j) = 5.7° provides large mass shift (~ 100 MeV) of the for the l^(j = 1/2) 
level and at the same time does not produce the mass shift of the 1^ level, which is almost pure 
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Table 2: The Ds(l^), Ds{2'^) meson mass shifts and widths due to the D*K decay channel for 
the mixing angle 4° (all in MeV) 



state 




^(theor) 


y^(cxp) 


T-,(theor) 
^ (D*K) 


-n ( eXD ) 

^ (D*K) 


5m 




2568(15) 


2458(15) 


2460 


X 


X 


-110 


Dsill) 


2537 


2535(15) 


2535(1) 


1.1 


< 1.3 


-2 




2575 


2573 


2573(2) 


0.03 


not seen 


-2 



Table 3: The 5^(1^)) -6^(2+) meson mass shifts and widths due to the B*K decay channel for 
the mixing angle 4° (all in MeV) 



state 


m(0) 


^(theor) 


^(exp) 


p(theor) 
^ {B*K) 


-p(exp) 
^ (B'K) 


5m 




5835(15) 


5727 


not seen 


X 


X 


-108 




5830(fit) 


5828 


5829 (1) 


0.8 


< 2.3 


-2 




5840(fit) 


5838 


5839(1) 


< 10"^ 


not seen 


-2 



j = I state. For illustration we show the scheme of the 1"*", 2+ shifts on Fig. [H We would like 
to stress here that the mass shifts weakly differ for Dg and Bg^ or weakly depend on the heavy 
quark mass: this can be directly illustrated using in the Eq. (|39p the expansion via the inverse 
heavy quark mass. 

6 Conclusions 

We have studied the mass shifts of the L>s(0+,1+ ) and 5^(0+, 1+ ) mesons due to strong 
coupling to the decay channels DK,D*K and BK,B*K. To this end the chiral quark-pion 
Lagrangian without fitting parameters has been used. 

We have shown that the emission of a NG meson, accompanied with the 75 factor, gives 
rise to maximal overlapping between the higher component with j = ^ of the P-wave meson 
{Ds,Bs) bispinor w.f. and the lower component (also with j = ^) of the S'-wave HL meson w.f. 
in considered 5- wave decay channel. Due to this effect, while taking the w.f. of the IP and 15" 
states with the use of the Dirac equation, large mass shifts of the 0+, 1+ states are obtained. 
In particular, the shifted masses M(S^,0+) = 5710(15) MeV and M{Bs, 1+') = 5730(15) MeV 
were calculated in agreement with the predictions in [H] and of S.Narison and by ~ 100 



Figure 4: Schemes of -05(1"*", 2"'") and 5^(1"^, 2"^) shifts due to chiral coupling 



B.dS 
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MeV lower than in [3], [3], [TO]. 

The widths of Dsi(2536) and i?si(5830) are also calculated. To satisfy the experimental 
condition r(Dsi(2536)) < 2.3 MeV the following limit on the mixing angle (j) (between the 
|j = I > and |j = ^ > states) is obtained: |(^| < 6°. 
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